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Summary. - The squared Laplace operator acting on symmetric rank-two tensor fields 
is studied on a (flat) Riemannian manifold with smooth boundary. Symmetry of this 
fourth-order elliptic operator is obtained provided that such tensor fields and their first (or 
second) normal derivatives are set to zero at the boundary. Strong ellipticity of the resulting 
boundary- value problems is also proved. Mixed boundary conditions are eventually studied 
which involve complementary projectors and tangential differential operators. In such a 
case, strong ellipticity is guaranteed if a pair of matrices are non-degenerate. These results 
find application to the analysis of quantum field theories on manifolds with boundary. 

PACS 03.70 - Theory of quantized fields. 

PACS 04.60 - Quantum gravity. 
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1. - Introduction 

Mathematicians and physicists have become familiar, along the years, with the key role 
played by the operators of Dirac and Laplace type in the investigation of elliptic operators 
on Riemannian manifolds. For example, it is by now well known that the symbol of the 
Dirac operator is a generator of all elliptic symbols on (closed) Riemannian manifolds 
[1], and deep results in index theory have been found by looking at non-local boundary 
conditions for operators of Dirac type [1,2]. More recently, local supersymmetry and 
quantum supergravity have led to the consideration of local boundary conditions for the 
Dirac operator [1,3-5], while theoretical models relevant for quantum chromo dynamics 
rely on mathematical structures weaker than the standard spin-structures where, again, 
a Dirac operator is found to be quite essential [1,6]. This analysis leads, in turn, to a 
deeper understanding of the geometry and topology of four-manifolds [6-8] . The operators 
of Laplace type, on the other hand, occur naturally in the consideration of quantized 
gauge theories (Abelian and non-Abelian) and Euclidean quantum gravity [5,9], and many 
efforts are devoted, within that framework, to the evaluation of the one-loop semiclassical 
approximation in quantum field theory, with the help of heat-kernel and ^-function methods 
[4,5,10,11]. 

Nevertheless, there is still room left for the analysis of many other classes of differential 
operators on manifolds. In particular, we are here concerned with the so-called conformally 
covariant operators [12,13]. They arise in the course of studying the behaviour of field 
theories and differential geometric objects under conformal rescalings of the background 
metric g. More precisely, if the conformal rescaling of g is written in the form g u = e 2aj g, 
a conformally covariant operator Q satisfies, by definition, the transformation property 

Q w = e -(-+4W2 Q( w = o) e (m-4) W / 2j (L1) 

where m is the dimension of the background Riemannian geometry. On compact Rieman- 
nian manifolds without boundary, the recent investigations have focused on the structure 
of anomalies [14] and on the application to the effective-action formalism [12,13]. For 
yet other classes of higher order differential operators, obtained by composition of oper- 
ators of Laplace type, much insight has been gained into the structure of the associated 
Green functions [15,16]. Even more recently, however, fourth-order elliptic operators have 
been studied on Riemannian four-manifolds with boundary, motivated by the path-integral 
quantization programme in gauges which are invariant under conformal rescalings [1,17- 
19]. Of course, even independently of the particular physical motivation, once that such 
operators have been studied for a first time in the presence of boundaries, the analysis 
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of a broader class of examples appears both natural and desirable, to improve the under- 
standing of some key features never studied before. For this purpose, bearing in mind 
all motivations mentioned so far, but leaving aside, for the time being, the immediate 
physical applications, we here consider the squared Laplace operator on flat Riemannian 
four-manifolds with smooth boundary, with its action on symmetric rank-two tensor fields 
on (M, g): 

U 2 = 9 ab 9 cd ^aV h V c Vd, (1.2) 
□ 2 : C°°(V,M) -> C°°(V,M). (1.3) 

With our notation, V is the connection on the vector bundle, V, of symmetric rank-two 
tensor fields on M. On the bundle V, we consider the particular DeWitt super-metric 
given by 



j^ab cd 



l(g ac g bd + g ad g bc )+ag ab g cd , (1.4) 



with inverse 



_i 1 / \ a 

E ab cd = g \9a C gbd + gadgbc ) - ( 1 + am \ gabg C d- (1.5) 



Since M is taken to be m-dimensional, with m^2, our form of the DeWitt super- metric 
is non-singular, and is positive-definite for all a > — — . 

The step leading to the definition (1.2) is simple but non-trivial, and hence deserves 
further comments. Indeed, in Euclidean quantum gravity at one loop, it is necessary to add 
to the Einstein-Hilbert action a gauge-averaging term, to ensure that the resulting gauge- 
field operator P ab cd admits a Green function and has a non-degenerate leading symbol 
[5,9]. In particular, the de Donder gauge-averaging term makes it possible to turn P ah cd 
into an operator of Laplace type (also called minimal), and in flat space one finds 

P ab cd h cd = -\Jh ab , (1.6) 

where 

□ =g cd V c V d = V e V e . (1.7) 

Thus, the proof of symmetry of P ab cd is reduced to proving that the operator Q defined in 
(1.7) is symmetric. In a flat Riemannian background, an analysis along similar lines leads 
to the consideration of the operator (1.2) starting from the composition of the operator 
P ab cd with itself, i.e. 

Pab Cd P cd ef hef=B 2 h ab , (1.8) 
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and hence the object of interest is actually the Q 2 operator (although we are leaving 
aside the problem of deriving the full gauge-field operator on manifolds with boundary for 
curvature-squared theories of gravity). 

Section 2 studies the inner product on the space of smooth and symmetric rank- 
two tensor fields on M, and performs the integrations by parts which are necessary to 
understand whether can be realized as a symmetric operator with suitable boundary 
conditions. A first set of boundary conditions, which make it possible to achieve the 
desired symmetry, are obtained in sect. 3. A second set of admissible boundary conditions 
in instead derived in sect. 4. Section 5 proves that the resulting boundary- value problems 
are strongly elliptic, and extends such a proof to the case of mixed boundary conditions 
for the squared Laplacian. Concluding remarks are presented in sect. 6. 

2. - Inner product and integration by parts 

In our paper, we use the following definition of inner product: 



where E ab cd is the DeWitt super- metric (1.4), and r) a },,hab are any two sections of the 
vector bundle of smooth, symmetric rank-two tensor fields over the m-dimensional Rie- 
mannian manifold (M,g). The definition (2.1) gives rise to an inner product if a+ ^ > 0, 
because by virtue of (1.4) one finds 



and, after using the general decomposition of rj a b, the coefficient of \J] a a J i n the integrand 
(2.2) is found to be a + — . Hereafter we shall therefore restrict a to be greater than — — . 



If Q 2 were a symmetric operator, the scalar products (jj, 2 /iJ and \ \^\ 2 rj,hj 

should be equal for all rj, h in the domain of Q 2 . In our paper, the above difference reads 
(hereafter, we assume M is four-dimensional) 




(2.1) 




(2.2) 
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- h ab V c V c V d V d {rj ab + ag ab r^j y/g d 4 x 

= J nP r]ooV p □ (/i 00 + ah) - h 00 V p □ (?7 00 + arj) 
+ 277oiV P D/i 0i -2h 0i V p nv° i 

+ VijV p \J (h ij + agVh) - h i:i V p \J {r] ij + ^7 d 3 x, 



(2.3) 



where we have assumed that g 00 = 1, g° l = 0, n p = (1,0,0,0) is the unit normal to the 
boundary, and 7 is the determinant of the induced three-metric on dM. Of course, h and 
rj in the integrand (2.3) denote the traces h = g ab h ab , r\ = g ab rj ab . It is now convenient to 
define, for <p ab = r) ab ,h ab: 



F 00 ( ip ) = ¥ > 00 + ag 00 <p = <p 00 + av, 

F 0i (i P ) = <P° i + <xg 0i V = <p \ 
F ij (ip) = (p ij +ag ij tp, 
so that Eq. (2.3) is expressed in the more convenient form 

Voo\JV P F 00 (h) - hoo\J^ P F 00 (v)] y/l d 3 x 
V0i\JV P F 0i (h) - h 0l \J^pF° l (v)] y/l d 3 x 
+ [ \ Vij \J^ P F ij (h) - □ Vpi^'fa)! d 3 x, 

JdM 1 J 



(2.4) 

(2.5) 
(2.6) 



JdM 



+ 2 



dM 



(2.7) 



where we have used the well known commutation property of covariant derivatives in a 
background whose Riemann curvature vanishes. Our aim is now to use integration by 
parts so as to remove the third-order derivatives occurring in Eq. (2.7). For this purpose, 
we use repeatedly the Leibniz rule. The analysis of the first boundary integral in Eq. (2.7) 
suggests considering (hereafter, Kij is the extrinsic-curvature tensor of the boundary, with 
trace K) 



(2.8) 



V p (nV Vm UF °{h)) = K V o o nF 00 (h)+nP(V pV oo)nF 00 (h) 

+ nV VQ0 \JV p F 00 (h), 
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jointly with 

V c (nP(V p ?7oo)V c F 00 (/i)) = K/(V p r ?00 )V c F 00 (/i) 
+ n p (V c V pVoo )V c F 00 (h) + n p (V pVoo )\jF 00 (h). (2.9) 

Moreover, the analysis of the second boundary integral in Eq. (2.7) suggests using the 
identity 

V p (n p VOl \jF 0l (h)) = K m \jF 0i (h) +n p (V pVoi )\jF 0i (h) 

+ n p mi \JV P F 0l (h), (2.10) 

and 

V c (n p (V p r] 0i )V c F 0i (hj) = K c P(V pmi )V c F 0i (h) 
+ n p (V c V p f]o l )V c F m (h)+n p (V pmi )nF°Hh). (2.11) 

Furthermore, the consideration of the third boundary integral in Eq. (2.7) makes it con- 
venient to use the identities 

V p {n p mj UF l3 {h)) = K Vij \jF ij (h) + n p (V pVij )\jF ij (h) 

+ n p VlJ \J^ P F lJ (h), (2.12) 



V c {n p (V pVij )V c FV(h)) = K p (V pVlJ )V c F^(h) 
+ n p (V c V pVlJ )V c F^(h)+n p (V pVlJ )\jF lJ (h). (2.13) 

If the boundary of M is smooth, one has ddM = 0. Moreover, the integral vanishes over 
zero-measure sets. It is hence possible to use Eqs. (2.8)-(2.13) and the Stokes theorem to 
re-express Eq. (2.7) in the form 



,B 2 h)-(B\h) 

-K^ooD^W-^ooDF 00 ^)) 
+ K? ((V^oo) V*F 00 (h) - (V J / i0 o)V J F 00 (r ? )) 
+ n p ((V c V p r ?00 )V c F 00 (/i) - (V C V p h 00 )V c F 00 ( V )^ 



-2K( VOl \jF 0t (h)-h 0l \jF 0l (v)) 

+ 2K l ^(V J ^ l F 0l (h) - (V^oi)V'F 0i (7/)) 

+ 2n p ([V c V pm )V c F 0i {h) - (V c V p /i 0l )V c F 0l (?7)) 

-K[ VlJ UF lJ (h)-h tJ UF lJ (v)) 

+ K{[{V r r lij )VF i ^h) - (V r /^)V^'(r/)) 

+ n p ((V c V p 7 /lJ )V c F lJ (/i) - (V C V P ^)V C F«(7/)) ^7 d 3 x. 



(2.14) 



3. - Boundary conditions: first option 

We know from the work of ref. [19], which studied the action of the squared Laplace 
operator on scalar functions / G C°°(M), that one can impose boundary conditions where 
/ and its normal derivative vanish at the boundary. This "doubling" of the boundary 
conditions, with respect to the analysis of the Laplacian, is clearly understood if one 
thinks of the eigenvalue problem. In other words, given a spectral resolution of a fourth- 
order elliptic operator, one deals with fourth-order eigenvalue equations which admit four 
linearly independent integrals. If one were to fix just the eigenfunctions at the boundary, 
one would not get enough equations to determine the coefficients of linear combination in 
the equation 



for the eigenfunction belonging to the eigenvalue A. At a deeper level, as shown in ref. [19], 
one has to integrate by parts to prove that suitable boundary conditions exist for which 
the operator is (essentially) self-adjoint (see Eq. (6.1)). 

In our problem, the technical details are more elaborated, since we study the action 
of Q 2 on smooth, symmetric rank- two tensor fields on a flat Riemannian manifold with 
smooth boundary (e.g. the Euclidean four-ball), but Eq. (2.14) can be used to derive 
all admissible sets of boundary conditions. First, we consider a scheme where half of the 
boundary conditions consist of requiring that all components of r\ a b and h a b, both spatial 
and normal, should vanish at the boundary. Thus, we require that 



fx(x) = Ci,a/i,aO*0 + C 2 ,xf2,x(x) + C 3 , x f 3 , x (x) + C 4 ,xh, x (x) 



Voo 



J dM 



J dM 



= 0, 



(3.1) 
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Voi 


dM 


hoi 


Vij 


dM 


h- 



dM 



dM 



= o, 



0. 



(3.2) 
(3.3) 



By virtue of (3.1)-(3.3) one then finds that the following covariant derivatives vanish at 
the boundary of M (see (3.20) and (3.21)): 



Vj?7oo 




VjTioo 




dM 








Vjh 0i 




dM 






dM 





dM 



= o, 



dM 



o, 



= 0. 



dM 



(3.4) 

(3.5) 
(3.6) 



At this stage, only the third, sixth and ninth line give non-vanishing contributions to the 
integrand in Eq. (2.14). It is hence appropriate to write them explicitly, bearing in mind, 
from sect. 2, that the normal to the boundary takes the form n p = (1,0,0,0). In other 
words, one has (with c,p ranging from through 3, and i, k ranging from 1 through 3) 

n p (V c V p r /00 )V c J P 00 (/i) = (VoV r ? oo)V F 00 (/i) + (V fc V r?oo) V k F 00 (h), (3.7) 

n p (V c V p r ?0l )V c F 0l (/i) = (V V ?7(k) V°F 0l (h) + (V fc V r ?0l )V fc J P 0l (/i), (3.8) 

n p {V c V pVij )V c F^{h) = (V V 7^)V i^(/i) + (V fc V ^)V fc F^'(/i), (3.9) 

and another triple of identities, which also contribute to the third, sixth and ninth line 
of (2.14), and are obtained by interchanging the roles of r\ a h and h a b in (3.7)-(3.9). The 
right-hand sides of (3.7)-(3.9) are linear combinations of products of covariant derivatives 
of rjab and h a b. Thus, it is sufficient to set to zero at the boundary only one of the two 
functions occurring in the product. For example, to avoid setting to zero on dM the second 
normal derivatives of r/ 0Q , r] 0i and rjij (and hence of /ioch ^ch an d hij as well) one can require 
that 

0, 



Voryoo 
V°F 00 (h)} = 

IdM 



dM 



dM 



dM 



dM 



Vo/lOO 

V /*oi 
'Vohij 
V°F 00 ( V ) 



dM 



dM 



= o, 



o, 



dM 



= o, 



(3.10) 
(3.11) 
(3.12) 
(3.13) 
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V°F 0i (h) 


dM 


V°F 0l (77)" 


dM 


= o, 


V°F lJ (h) 


dM 




dM 


= 0. 



(3.14) 
(3.15) 



Note that we do not get an overdetermined problem, because, by virtue of (2.4)-(2.6) and 
(3.10)-(3.12), Eqs. (3.13)-(3.15) are satisfied. Now it is helpful to consider an example, 
and for this purpose we choose the Euclidean four-ball, whose metric may be locally cast 
in the form 



g = dr <E> dr + r 2 Cijdx l ® dx\ 



(3.16) 



where the radial coordinate r lies in the closed interval [0,a], with a the radius of the 
three-sphere boundary, and Cy is the metric on a unit three-sphere, with local coordinates 
{x % }. One then finds (Latin indices run here from 1 through 3) 



Vo/ioo — 



dh 



oo 



dr ' 



V7 U dfl0i 1 h 
Vo^Oi — n 0i, 

Or t 
dhij 2 

V7 h dh °° 2 h 

Vfcftoo — z ^Ofc, 

OHq' 1 1 
Vfc/ioi = -rnr - r ki hi h ik H — g ik h 



00 1 



dx k T T 

with T used to denote the connection coefficients, and hence 



^ 9 9/ioo 1 /<9/i o . O <9^ofc 4 
VfcVo/ioo = ^— + z ~a n 0k 

ox K or r \ ax K or r 



(3.17) 

(3.18) 
(3.19) 
(3.20) 
(3.21) 



(3.22) 



VfcV /ioi 



9 / dh i 1 



h 



Oi 



1 

T 



0i 



r i (dhoi 1 VI 



dx k 

dh Q0 
dr 



^00 

r 



1 

r 



dr 



-hit 
r 



(3.23) 
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V fc Vo/i 



d 

dx k 



-^r 1 ~ -hij) - -Vfc/j, 
or t I r 



2 ( d 1, 



i)o 



(3.24) 



The joint effect of Eqs. (3.1)-(3.3), (3.7)-(3.12) is then that the right-hand side of Eq. 
(2.14) vanishes if the following boundary conditions hold: 



h 



ab 



dM 



Vab 



dM 



= o Va, 6= 0,1,2,3, 



~dh ab ~ 




dt]ab 




. dr _ 


dM 


. dr _ 


dM 



Va, b = 0,1, 2, 3. 



(3.25) 
(3.26a) 



One then deals with 10 boundary conditions on h a b,r} a b and 10 boundary conditions on 
their normal derivatives, bearing in mind the symmetry of these rank-two tensor fields. 
This is a generalization of the boundary conditions obtained in ref. [19] for the squared 
Laplacian acting on smooth functions, and once more the number of boundary conditions 
is exactly doubled, with respect to the analysis of the Laplacian. Last, a covariant form of 
Eq. (3.26a) is easily obtained from (3.10)-(3.12), i.e. 



n p Vph ab 



J dM 



n p VM = Va, b = 0,1,2,3. 

' ' dM 



(3.266) 



The normal components of <p a b = h a b or r\ a b are given by 



(3.27) 



V?oo = VabU n 

Voi = ¥ain a , (3.28) 
and the spatial components are obtained, instead, by applying a projection operator, i.e. 



where, on defining [20] 
one has [20] 



q a b = 5 a b - n a n\ 



ncd „ c „ d 
„h = 



(3.29) 

(3.30) 
(3.31) 



l ab = Q(a %) 

The tensor field q b is the standard projector of tensor fields over the bounding surface, 
and is, by construction, orthogonal to the unit normal vector, in that q b nb = n a — n a = 0. 
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4. - Boundary conditions: second option 



After writing down the identities (3.7)-(3.9), one can also make the alternative choice, 
according to which (cf. (3.10)-(3.15)) 



VoVo^oo 

V V ??oi 
V V ??ij 
V k F 00 (h) 
V k F 0l (h) 
V k F ij (h) 



J dM 



dM 



dM 



dM 



J dM 



dM 



VoVo/lQO 



J dM 



VoVo/iOi 

'Vo'Vohij 
V k F 00 ( v ) 
V k F 0i ( v ) 
V k F ij (rj) 



dM 



dM 



dM 



dM 



dM 



= o, 

= 0, 
= 0, 
= 0, 
= 0, 
= 0. 



(4.1) 

(4.2) 
(4.3) 
(4.4) 
(4.5) 
(4.6) 



Now the boundary conditions (3.1)-(3.3), jointly with the definitions (2.4)-(2.6), imply 
that Eqs. (4.4)-(4.6) are identically satisfied. In particular, on the Euclidean four-ball, by 
virtue of the formulae 



VoVq/ioo 



d 2 h 



00 



VoVo/iOi — 



d 2 h 0i _ 2 dhpj 
dr 2 t dr 



H nh 0i , 



8 2 h- A8h- fi 



(4.7) 
(4.8) 
(4.9) 



the boundary conditions (3.1)-(3.3) and (4.1)-(4.3) lead to 10 homogeneous Dirichlet con- 
ditions, jointly with 10 boundary conditions for a linear combination of first- and second- 
order partial derivatives, i.e. 



(4.10) 

= 0, (4.11) 
0. (4.12) 





'd 2 r]oo' 
. dr 2 _ 


dM 


~d 2 h o~ 
. dr 2 _ 


= 0, 

dM 




~d 2 rioi 


2 dr/oi ' 




'd 2 h 0i 


2 dh 0i ' 




Or 2 


r dr 


dM 


_ dr 2 


r dr 


dM 


'd 2 f] l3 


4 druj ' 




~d 2 hij 


Adhif 




_ dr 2 


t dr 


dM 


_ 8t 2 


t dr 


dM 
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Of course, the boundary conditions (4.1)-(4.3) can also be expressed in a covariant form 
by a single equation, i.e. 



n p n q V p V q r] a i 



dM 



n p n q V p V„h ab 



J dM 



OVo, 6= 0,1,2,3. 



(4.13) 



As in sect. 3, the number of boundary conditions is doubled with respect to the analysis 
of the Laplacian, and Eqs. (3.1)-(3.3) and (4.1)-(4.3) provide the generalization of yet 
another set of boundary conditions derived, in ref. [19], when the action of the squared 
Laplacian on smooth functions is considered. 



5. - Strong ellipticity 

If a differential operator is studied on a Riemannian manifold M with smooth bound- 
ary dM, the ellipticity is obtained upon proving that the leading symbol is elliptic in the 
interior of M, and that a unique solution exists of the eigenvalue equation for the leading 
symbol, subject to a decay condition at infinity and to suitable boundary conditions. For 
the squared Laplace operator defined in Eq. (1.2), the leading symbol 

^(□ 2 ;*,£) = l£l 4 ^ = & 

with £ G T*(M), is clearly elliptic. In fact, for any Riemannian metric g, the leading 
symbol in (5.1) is non-degenerate for £ ^ 0. Moreover, for any A G C — 3?+, one has 

/ / \ \ / \ dim(V) 

det(a L (D 2 ;x,e) -Aj = (|£| 4 - A] ? 0. (5.2) 

The left-hand side of Eq. (5.2) vanishes if and only if both £ and A vanish. This completes 
the proof of ellipticity of cl(0 2 ) i n the interior of M. 

Now we have to define the condition of strong ellipticity, which involves in a crucial 
way the boundary conditions. For this purpose, we consider again the leading symbol 

<jl ^O 2 ; x, r, (, uj, with |x fc } local coordinates on dM, r the normal geodesic distance 

to the boundary, Q as components of a cotangent vector on the boundary, and u a real 
parameter. We then set r = and replace to by — id r , and consider the eigenvalue equation 
for a L , i.e. [9,10] 

a L (D 2 ; x, r = 0, C, -id r )<f(r) = \<p(r). (5.3) 
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Note that, in writing down Eq. (5.3), we have not doubled the number of arguments in 
the leading symbol. In the presence of boundaries, the m local coordinates x are split into 
m — 1 local coordinates on the boundary, jointly with the normal geodesic distance (which, 
of course, vanishes on the boundary). Similarly, the m coordinates £^ are split into m — 1 
coordinates on T*(<9M), appropriate for cotangent vectors, jointly with a real parameter 
uj. Hence one has (with |C| 2 = CiC) 

^(□ 2 ; {x k } , r = 0, {0} ,u) = |C| 4 + ^ + 2u 2 \(\ 2 . 

Last, to obtain the left-hand side of Eq. (5.3), one has to replace u by — id r as we said 
before [10], and this leads to the fourth-order differential operator 

^(□ 2 ; {x k } , r = 0, {Q} , co = -id r ) = ^ - 2|C| 2 ^ + |C| 4 - 

Thus, Eq. (5.3) should be viewed as an ordinary differential equation involving the field 
if which is a smooth section of the vector bundle V over the product manifold dM x 3? + . 
By definition, for a given form of the boundary operator B, the boundary- value problem 

^ Q 2 , B j is strongly elliptic with respect to the cone C — K + if there exists a unique 

solution of Eq. (5.3) with A G C — 3? + , for all (£, A) 7^ (0,0), subject to the asymptotic 
condition 

lim (p(r) = 0, (5.4) 

1 — >oo 

and to the boundary condition 

a g (B)(x,(Mv) = i>', (5.5) 

for all ip'. With a standard notation, ip(<p) denotes the boundary data, usually arranged in 
the form of a column vector consisting of the field and its normal derivative(s) evaluated 
at the boundary. The boundary conditions originally imposed read (cf. (3.25) and (3.26), 
or (3.25) and (4.13)) 

Bip(<p) = at dM, (5.6) 
where, in our case, the boundary operator takes the form 
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Unlike the case of mixed boundary conditions for the Laplacian, no projectors occur in 
Eq. (5.7), because we are setting to zero on the boundary the whole of ip and its first (or 
second) normal derivative. The graded leading symbol cr g (B) of B, is defined according to 
the rule [10] 

cr g {B) i:j = a L (B) i:j if ord(S y ) = j - i, (5.8a) 

<r g {B)ij = if ord(Sy) <j-i, (5.86) 

where "ord" denotes the order of the Bij component of B as a differential operator on the 
boundary, subject to the restriction 

ord(By) <j-i. (5.9) 

By virtue of (5.7)-(5.9), the graded leading symbol of the boundary operator is again the 
identity matrix 

°9(B)=( I f) ( 5 - 10 ) 

in the boundary-value problems of Sees. Ill and IV, whereas ip((f) consists of the column 
vectors (here, the symbol -n coincides with Vo used in the previous sections) 

^)=( Mom 



[<P; 



N\dM 



or 

Vp]dM 



V VP\NN\dM 

respectively, with ip' any pair of boundary data [9] 



After having defined in detail the strong ellipticity setting, we can perform the next 
step, i.e. the solution of Eq. (5.3) jointly with the conditions (5.4) and (5.5). Indeed, Eq. 
(5.3) reads 



(p(r) = X<p(r). (5.11) 
Equation (5.11) is solved by (p(r) = e ar , with a given by the roots of the algebraic equation 

a 4 -2a 2 |C| 2 + |C| 4 -A = 0, (5.12) 
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i.e. 



«i = +y / KI 2 + v / A, 

« 2 = +V / KI 2 -v / A, 
«3 = -V / KI 2 + v / A, 
«4 = -V / KI 2 -v / A. 



(5.13) 
(5.14) 
(5.15) 
(5.16) 



Comparison with the case of the Laplace operator [10] shows that one obtains strong 
ellipticity provided that ±\/A G C — 9ft+, which yields 



AG (C-3?+)nC = C-K+. 



(5.17) 



Among the four values of a, only «3 and «4 fulfill the condition (5.4). The exponentially 
decaying solution picked out by (5.4) reads therefore (xi an d X2 being some parameters) 



where 



One then finds 



(p(r) = xie pir + X2e 



Pl = +V / KI 2 + v / A, 
P2 = +V / KI 2 -v / A. 



dtp 
dr 



(p\r=0 = Xl +X2, 

= -PiXi - P2X21 



r=0 



d 2 if 



dr 2 



= P1X1 + P2X2- 



(5.18) 

(5.19) 
(5.20) 

(5.21) 
(5.22) 

(5.23) 



r=0 



Thus, if the boundary conditions (3.25) and (3.26) are imposed, the condition (5.5) leads 
to 

I 0\ / [(p]dM \ = (ip'o 

1 J \ MdMj W[ 



(5.24) 



which implies, by virtue of (5.21) and (5.22), 



Ax = 
15 



(5.25) 



where y is the column vector ( ^ 1 | , and A is the matrix 

\X2j 

A= ( 1 1 V (5.26) 

V -pi -pi) 

Now a unique solution x = A~ lr ip' exists of the linear system (5.25), because 

det(A) = Pl -p 2 ^0. (5.27) 

Explicitly, one finds 

X. - , (5.28) 

(p2 - Pi) 

X2 = JA±flf, (5 . 29) 

(P2 - Pi) 

Moreover, if the boundary conditions (3.25) and (4.13) are imposed, the condition (5.5) 
leads to 

I o\ / [<p] eM \ = U' Q 

I J \ \}P\N n\&m J 



(5.30) 



A = ( \ \ ) • (5-32) 



and this implies, by virtue of (5.21) and (5.23), 

A X = tf, (5-31) 

where A is the matrix 

.Pi P2 

A unique solution x = A~ l ip' exists of the linear system (5.31), because 

det(A) = p\ - p\ = -2V\ ^ 0. (5.33) 
The explicit calculation yields (cf. (5.28) and (5.29)) 

_ (^ - pMI (5 M) 

X1 ~ (Pi -Pi) ' (5 - 3 ) 

X2 = - iPl-Pl) • (5 - 35) 



16 



So far our examples of strong ellipticity have been almost straightforward. A more 
relevant case is obtained on considering mixed boundary conditions for the squared Lapla- 
cian, with boundary operator involving also tangential derivatives. For this purpose, we 
assume that the boundary operator in (5.7) is replaced by 

B=(? k Al (5-36) 



N A fc I-IL, 

where k is an integer > 1, A is a first-order tangential differential operator on the boundary: 

A : C°°(W , dM) -> C°°(W , dM), 

and IT is a self-adjoint projector. With a standard notation, W is the bundle of boundary 
data over dM, given by the direct sum 

W = W ®W U (5.37) 

with Wj representing normal derivatives of order j. In the previous case, we had W 
expressed by (5.37) or, instead, by 

W = W Q ®W 2 . (5.38) 

The projector IT maps each Wk sub-space into itself. The form of A which leads to tan- 
gential derivatives in the boundary conditions is 



A=(/-n) 



^(rVi + Vir*) +s 



(J -II), (5.39) 



where T l are endomorphism- valued vector fields on dM, and S is an endomorphism of 
the vector bundle W . Following ref. 9, T l and S are taken to be anti-self-adjoint and 
self-adjoint, respectively, and annihilated by the projector IT from the left and from the 
right. By virtue of the assumption (5.36), the graded leading symbol of the boundary 
operator is (here T = T^Q) 

= ((")' /-n)- < 5 - 40 > 

and hence, when (5.37) holds, the boundary condition (5.5) takes the form 

n _ \( X1+X2 (541) 



(iT) k I-Uj \-piXi-P2X2 UW' 
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where 

$> = n(7i + 72), (5-42) 
^ = (I-n)(-p l7l -p 2 7 2 ). (5.43) 

This means that the right-hand side of Eq. (5.41) is a smooth section of an auxiliary vector 
bundle W over <9M, where any ip' e C°°(W, 9M) is given by ip' = Pip, with 

,° n ). (5.44, 

and -0 G C°°(W, dM). The boundary operator in (5.36) is related to the operator P defined 
in (5.44) by the equation 

B = PL, (5.45) 



where 

' I 
A fc I 



(5.46) 



In other words, since we deal with a squared Laplacian, we allow for a number of tangential 
derivatives greater than the number considered in the case of the Laplacian, and we try 
to impose mixed boundary conditions, so that projectors (rather than identity operators) 
occur in the boundary operator. Since ip differs in general from ip, we have constants 71 
and 72 in (5.42) and (5.43) instead of xi and X2 (cf. (5.18)). 

Thus, for all j = 1, 2, if the matrix (iT) k —Ipj is non-singular we may use the identity 

Xj =U Xj + (I-U) Xj (5.47) 

to solve the system (5.41) in the form 

n Xj = n 7i , (5.48) 

(I - U) Xj =(l Pj - {iT) k ) _1 [(zT) k U + (I - U) Pj ] 7j , (5.49) 

A necessary and sufficient condition to obtain a unique solution is of course the non- 
degeneracy of the matrix (Ipj — (iT) k ) : i.e. 



det 



I Pj -(iT) k ^0. (5.50) 
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Note now that (see (5.19) and (5.20)) 



I Pl -(iT) k = IyJ\C\ 2 + V\-(iT) k , 



(5.51) 



iP2-m k = iy/\c\ 2 ->ft-m k , 



(5.52) 



and such matrices are never singular if the condition (5.17) is satisfied and iT is self-adjoint. 

6. - Concluding remarks 

Our paper has analyzed the squared Laplace operator Q 2 acting on symmetric rank- 
two tensor fields ip a b on (flat) Riemannian manifolds with boundary. Its original contribu- 
tions, of technical nature, are as follows. 

(i) Symmetry of Q 2 is achieved provided that both Lp a b and its normal derivative n p Vpip a b, 
or ip a b and the second normal derivative n p n q V P V q <~p a b, are set to zero at the boundary. 

(ii) The resulting boundary-value problems are strongly elliptic with respect to the cone 



(iii) Strong ellipticity with respect to C — 3?+, in the case of mixed boundary conditions in- 
cluding tangential derivatives, has also been proved. Interestingly, no restriction involving 
T, and hence the matrices T- 7 , is obtained unlike the case of an operator of Laplace type 



Of course, in the case of symmetric rank-two tensor fields the identity operator in Eq. 
(5.25) reads actually 5^ a <5&, but apart from such minor details, all calculations in Sec. V 

prove indeed that the boundary- value problem ^0 2 ,-B^ is strongly elliptic with respect 

to the cone C — 3?+, on considering the vector bundle of symmetric rank-two tensor fields 
over M. Our proof is simple but of some interest, because it clearly shows the role played 
by fourth-order operators in doubling the number of linearly independent solutions of the 
eigenvalue equation (5.3) for the leading symbol. Strong ellipticity is crucial to ensure the 
existence of the asymptotic expansions used in the theory of heat-kernel asymptotics [10]. 
From the point of view of quantum field theory, this means that the one-loop semiclassical 
approximation is well defined and can be explicitly evaluated [19]. 

We find it appropriate to stress once more that the double integration by parts used 
to derive Eq. (2.14) is necessary to recover the correct number (and form) of boundary 
conditions for a fourth-order elliptic operator like Q 2 . For example, its simplest (but still 



[9]. 
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,4 

useful) form, i.e. the operator B = on a closed interval of the real line, satisfies the 
identity [19] 



(Bu,v) - {u,B ] v) 



d 3 



u 



dx 3 



d z u* dv 
dx 2 dx 



-, 1 



+ 



du* d 2 v 
dx dx 2 



u 



dx 3 



n 1 



(6.1) 



where u is a vector in the domain of B, v is a vector in the domain of the adjoint B^ of 
B, and we use the definition of inner product [19] 



(u,v) 



u*(x)v(x)dx. 



(6.2) 



Domains of self adjointness of are therefore, in particular, the set of functions belonging 
to AC 4 [0, 1] (this is the set of functions in L 2 [0, 1] whose weak derivatives up to third order 
are absolutely continuous in [0,1], which ensures that the weak derivatives, up to fourth 
order, are Lebesgue summable in [0,1], and that all u in the domain are of class C 4 on 
[0,1]) and satisfying the boundary conditions [19] (cf. (3.25) and (3.26)) 



u(0) = u(l) = 0, 

u'(0) = u'(l) = 0, 
or Eq. (6.3) jointly with [19] (cf. (3.25) and (4.13)) 

w"(0) = u"{\) = 0. 



(6.3) 
(6.4) 

(6.5) 



Thus, a complete correspondence can be established between the boundary-value problem 
for the squared Laplace operator in one dimension [19] and the more elaborated case 
studied in our paper. 

It now appears both interesting and necessary to study a scheme more general than 
the one where 7] ab ,h ab and their first or second normal derivatives are set to zero at the 
boundary. For this purpose, one may start again from Eq. (2.14), but requiring that the 
integrand, as a whole, should vanish at dM. Such an integrand a(r), h) may be re-expressed 



as 



a( V , h)=[- Kh ab V c + nP{V c V p h ah ) + KP(W p h ab )\ V c F ab {r,) 

- K Vab V c + nP(V c V pVab ) + K c P(V pVab )] V c F ab (h). (6.6) 
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The non-trivial problem, however, is to understand how to derive the boundary opera- 
tor (5.36) from the analysis of (6.6). In the previous section, the introduction of (5.36) 
was too simplified, leaving aside the problem of integrating by parts in the action. If it 
were possible to achieve this, it would then remain to be seen whether such a kind of 
generalized boundary conditions for Q 2 can be derived from an invariance principle, as is 
indeed the case for the Laplacian itself [1,5,9], upon requiring invariance under infinitesi- 
mal gauge transformations. Non-local boundary conditions for Q 2 might also be studied 
with some profit, following the recent attempts to consider a non-local formulation of Eu- 
clidean quantum gravity based on integro-differential boundary conditions [21]. Last, but 
not least, the resulting heat-kernel asymptotics should be thoroughly developed, to supple- 
ment the recent, encouraging progress in the case of generalized boundary-value problems 
for operators of Laplace type [22-25]. All this adds evidence in favour of the problems 
of quantum field theory and spectral geometry being able to lead to a deeper vision in 
modern mathematical physics [1]. 

* * * 

The author is much indebted to Ivan Avramidi for scientific collaboration on the tech- 
niques used in the present paper and for comments which led to a substantial improvement. 
This work has been partially supported by PRIN97 "Sintesi" . 
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